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INTRODUCTION j 

Tlii3 final report summarizes the research accomplished under NASA 
Crant NCR 16-002-038 during the entire funded period which extended from 
November 15, 1971 to March 31, 198A. The total funding amounted to 
$261,160. The Technical Officers for this grant were the late John V. 

Rakich and Joe C. Marvin of NASA Ames Research Center. 

The research performed during the grant period can be divided into 
the following major areas: * 

I 

1. Development of equilibrium air curve fits t 

2. Computation of hypersonic rarefied leading edge flows 

3. Computation of 2-D and 3-D blunt body laminar flows with an 
impinging shock 

4. Development of a two-dimensional or axisymmetric real gas blunt 
body code 

5. Study of an over-relaxation procedure for the MacCormack finite- 
difference scheme 

6. Computation of 2-D blunt body turbulent flows wi.h an impinging 
shock 

7. Computation of supersonic viscous flow over delta wings at high 
angles of attack 

8. Development of a generalized PNS code 

9. Computation of the Space Shuttle Orbiter flowfield 

These major areas of research are discussed separately in the next section. 
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DESCRIPTION OF PERFORMED RESEARCH 

1. Development of equilibrium air curve fits 

The thermodynamic and transport properties of equilibrium air are 
required for the numerical computation of flowfields around bodies in 
high-speed flight. In particular, when using a finite-difference tech- 
nique to solve the unsteady Navi r-Stokes equations in conservative form, 
it becomes necessary to determine the thermodynamic properties as a 
function of density (0) and internal energy (e) . This requirement led 
to the present study in which two different approaches were devised for 
the case of equilibrium air. In the first approach, the NASA-Ames RGAS 
program was modified* to allow density and internal energy to be the 
independent variables. This approach permits a very accurate determina- 
tion of the thermodynamic properties of air. Unfortunately, the table- 
lookup feature of the RGAS program is too cumbersome to be effectively 
employed on advanced computers. For this reason, and also to reduce 

computation time on conventional serial computers, simplified curve fits 
2-3 

were developed in the second approach. 

The simplified curve fits were constructed using Grabau-type transi- 
tion functions. A transition function of this type can be used to 
smoothly connect two surfaces together with or without a line of inflec- 
tion between them. The coefficients In the transition functions were 
determined using a least squares computer program to fit the data from 
the NASA RGAS program. The resulting curve fits include the following 
correlations: 
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p - p(e,o) T - T(p,p) 

T - T(e,p) U “ U(e,p) 

a - a(e, p) k - k(c,p) 

h - h(p,p) 

The ranges of validity for these curve fits are the same as the NASA RGAS 

program, namely, temperatures up to 25,000°K and densities from 10 ^ to 
-3 

10 amagats (i.e., p/p q ). The maximum percentage differences between 
the simplified curve fits and the RCAS program for the primary variables 
p * p(e,t!) and h ■ h(p,p) are 4.7% and 4.6%, respectively. 


Computation of hypersonic rarefied leading edge flows 
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The hypersonic leading edge problem has been the subject of a con- 
siderable number of both experimental and theoretical investigations over 
the past three decades. The theoretical investigations can be divided 
into two groups, depending on whether a kinetic theory approach or a 
continuum flow approach is used. Prior to the present study, the con- 
tinuum approach was primarily based on solving simplified forms of the 
complete Kavier-Stokes equations. The accuracy of these solutions was 
somewhat questionable because of the neglected terms. In the present 

study, the hypersonic rarefied flow near the sharp leading edge of a 

4-6 

flat plate was computed using a finite-difference solution of the 
complete unsteady Navier-Stokes equations. MacCormack's original ex- 
plicit, second-order accurate, finite-difference scheme was used to solve 
the equations at each interior grid point. The computational region 
extended from the freestream ahead of the plate to the strong-interaction 
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regime. Both wall slip and temperature jump were included in the cal- 
culations. The numerical results compared favorably with experimental 
data and Monte Carlo simulations. In addition, it was shown for the 
first time that a numerical solution of the Navier-Stokes equations 
predicts that the wall pressure asymptotically approaches the free-mole- 
cule limit at the leading edge. 

7—8 

In addition to this study, a related investigation was conducted 
to determine if a numerical solution of the complete Burnett equations 
would yield more accurate results for the hypersonic leading edge problem. 
The Burnett equations are a higher-order set of equations (in terms of 
Knudsen number) than the Navier-Stokes equations, and thus it might be 
reasonable to assume that they would yield better results than the Navier- 
Stokes equations in rarefied regions. However, it is believed by many 
that while the Burnett equations can be used to improve the solution in 
regions where the Navier-Stokes equations give good accui«cy, it is not 
safe to assume that they can be used to progress into more rarefied 
regions where the Navier-Stokes equations are invalid. The present inves- 
tigation was undertaken to test this latter assertion. The complete 
Burnett equations were integrated using MacCormack’s explicit finite- 
difference scheme. The resulting solution was compared with experimental 
data, a numerical solution of the complete Navier-Stokes equations, and 
a Monte Carlo simulation. Based on this comparison, it was evident that 
the Burnett equations in conjunction with Schamberg's boundary condition 
gave a much less accurate description of the rarefied flowfield near the 
leading edge of the flat plate. This result tends to confirm the belief 
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that the Burnett equations cannot be used to progress into rarefied 
regions where the Navier-Stokes equations are already invalid. 

3. Computation of 2-D and 3-D blunt body laminar flows with an impinging 
shock 

An extraneous shock wave impinging on a blunt body in a hypersonic 
flow has been observed to greatly increase both the heat transfer rate 
and pressure near the impingement point. Flowfields of this type occur 
on the Space Shuttle and other maneuverable re-entry vehicles. The in- 
tense heating and high pressures occur over a small region where a dis- 
turbance. originating at the intersection of the impinging shock and bow 
shock, strikes the body. The disturbance may be a free shear layer, a 
supersonic jet, or a shock wave depending on the strength and location 
of the impinging shock and the shape of the body. Edney has described 
six different types of shock interference patterns which can occur. 

Because of the very complicated nature of shock impingement flow- 
fields, previous attempts at predicting the maximum heating rates and 
pressures were limited tj scmiempirlcal approaches. In the present study, 
this empiricism was eliminated by numerically computing the entire shock 
impingement flowfield. This was accomplished by using a "time-dependent," 
finite-difference method to solve the complete set of compressible Navier- 
Stokes equations for a laminar flow. The major reason for using the 
"time-dependent" method is that the resulting unsteady Navier-Stokes 
equations are a mixed set of hyperbolic-parabolic equations for both sub- 
sonic and supersonic flows. As a result, a very complicated shock 
impingement flowfield which contains both subsonic and supersonic regions 
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can be calculated as an initial-value problem. An additional advantage 
is that since the complete Navier-Stokes equations are solved, all shock 
waves, shear layers, jets, and the wall boundary layer are "captured" 
automatically in the solution without prior knowledge of their location 
or even existence. 

For very low Reynolds number flows, the entire shock impingement 

flowfield, including the bow shock, can be "captured" using this so-called 

"shock-capturing" approach. This approach was used for computations made 

9 

during the early stages of the present study. However, it was found that 
for higher Reynolds number flows it was not practical to "capture" the 
bow shock because of the numerical difficulties associated with the large 
gradients at the bow shock. Instead, it was found to be more convenient 
to treat the bow shock as a discontinuity, across which the Rankine- 
Hugoniot equations could be applied, while leaving the boundary layer and 
interaction regions to be captured as before. This latter approach is 
the so-called "shock-fitting" method. 

MacCormack's original, explicit, finite-difference algorithm was 
used to solve the Navier-Stokes equations at each interior grid point. 
Initially, the blunt body flow without the impinging shock was computed. 
The impinging shock was then introduced into the flowfield by specifying 
its location and strength. The computation was then restarted and con- 
tinued until the final "steady-state" solution was reached. 

Both two-dimensional and three-dimensional blunt body flows with 
planar impinging shocks were computed. In the two-dimensional cases, 

the planar impinging shock wa3 parallel to the centerline of the infinite 

9-11 

cylinder. Early 2-D results were qualitatively compared with existing 
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3-D experiments because of the lack of suitable 2-D experiments. Later. 

it became possible to compare 2-D computations with 2-D experiments con- 

12 

ducted by J. W. Keyes at NASA Langley. The numerical results were in 

excellent agreement with the 2-D experimental values of shock shape. 

13-15 

pressures, and heat transfer rates. The 3-D computations, showed 

reasonable agreement with the experimental results. These computations 
were hampered by the lack of grid points which was the direct result of 
computer storage limitations. 

4. Development of a two-dimensional or axisymmctric, real gas, blun t 
body code 

A two-dimensional or axisymmetric, viscous, real gas, blunt body code 
3 

was developed in this study. The complete Navier-Stokes equations in 
curvilinear, orthogonal coordinates were employed as the governing equa- 
tions. The use of orthogonal coordinates simplifies the application of 
body boundary conditions and the determination of a starting line for 
supersonic calculations downstream of the nose. The steady-state solu- 
tion to the governing equations is obtained using the original, explicit 
MacCormack scheme. The bow shock is treated as a discontinuity using 
the Rankine-Hugoniot relations. 

Both perfect gas and equilibrium air calculations can be performed 
with this code. The equilibrium air flow properties are introduced into 
the calculation procedure via the simplified curve' fits described earlier. 
This code was tested by calculating both perfect gas and equilibrium air 
flows over standard sphere test cases. 


5. Study of an over-relaxation procedure for the MacCormack finite- 
difference scheme 

An over-relaxation procedure was developed*^ ^ for the MacCormack 
finite-difference scheme in order to reduce the computation time required 
to obtain a steady-state solution. The implementation of this accelera- 
tion procedure to an existing computer program using the original explicit 
MacCormack method is extremely simple and does not require additional 
storage. The over-relaxation procedure does not alter the steady-state 
solution, which is second-order accurate. The method was first applied 
to Burgers' equation. A stability condition and an expression for the 
increase in the rate of convergence was derived. The method was then 
applied to the calculation of the hypersonic viscous flow over a flat 
plate using the complete Navier-Stokes equations, and the invlscid flow 
over a wedge using the Euler equations. Reductions in computing time by 
factors of 3 and 1.5, respectively, were obtained by the over-relaxation 
procedure. 

6. Computation of 2-D blunt body turbulent flows with an impinging 3hock 

During this study, a computer program for computing two-dimensional 
turbulent blunt body flows with an impinging shock wave was developed. 

This code utilizes the time-dependent, time-averaged Navier-Stokes equations 
with algebraic eddy viscosity and turbulent Prandtl number models emoloyed 
in the expressions for shear stress and heat flux. For the turbulent 
boundary layer, the two-layer mixing length model of Cebeci and Smith is 
used. This model contains a correction for the streamwise pressure 
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gradient which occurs in blunt body ilows. For the turbulent shear layer 
a modified Prandtl mixing length model is used. In the region where the 
shear layer impinges on the boundary layer, a combination of the two 
models is employed. 

Originally, the standard, unsplit, MacCormack finite-difference 
scheme was used to integrate the governing equations. Tills proved unsuc- 
cessful, however, because the allowable time step is proportional to the 
grid spacing. Since, a very fine mish is required to prcperly resolve 
the turbulent boundary layer it was found that excessively long computer 
times were necessary to obtain a steady-state solution. Because of this 
difficulty, the Beam-Warming implicit scheme was incorporated into the 
turbulent code. This scheme increases the computation time per step by 
a factor of 1.6 but permits a time step which is many times greater than 
the time step of the previous explicit scheme. Thus, steady-state 
turbulent solutions can be obtained with substantially fewer iterations. 
The present implicit scheme is noniterative, spatially factored, and 
employes Euler time differencing which makes the scheme second-order 
accurate in space but first-order accurate in time. For the present 
steady flows, time accuracy is not important. The irregular-shaped bow 
shock is treated as a discontinuity across which the Rankine-Hugoniot 
equations are applied. The code was used^ ^ to compute a Type III 
turbulent shock interference flowfield. The computed results were in 
excellent agreement with the experimental values of shock shape, wall 
pressures, and heat transfer rates. 



In a related study! the present code (with the MacOriuack finite- 

20 

difference scheme) was used to calculate the turbulent hypersonic flow 
over severely indented nose tips. The computed results were compared 
with experimental data from wind-tunr.el tests. Qualitative agreement was 
obtained for the surface pressure distribution and flowfield structure! 
including the separated bubble in the indented region. 

7 . Computation of supersonic viscous flow over delta wings at high 
a rgles of attack 

21-23 

Two complementary procedures were developed to calculate the 

supersonic, viscous flow over conical shapes a„ jarge angles of attack, 
with application to cones and delta wings. In the first approach, the 
flow was assumed to be conically self-similar. This approximatic i is 
suggested by the results of experiments for supersonic flows around coni- 
cal bodies and wings. The resulting Navier-Stokes equations are solved 
at a given Reynolds number with a time-marching explicit f inite-dif rcrence 
algorithm. These equations are referred to as rhe conical Navier-Stokes 
equations. These equations had been used previously to calculate the flow 
around cones at high angles of attack and delta wings with supersonic 
leading edges, however, in those calculations the bow shock was captured 
and the body shapes were limited by a rather restrictive geometry. The 
present method treats the shock as a sharp discontinuity and allows for 
a completely general crcss-sectional shape and distribution of the finite- 
difference grid points. 

In the second approach, only the streamwise viscous derivatives were 
neglected in the steady Navier-Stokes equations. The resulting coe- Cions 
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are referred to as the "parabolized" Navier-ftokes (PNS) equations. The 
PNS equations are hyperbolic-parabolic with respect to the streamwlse 
direction. This permits the solution ro be marched downstream from a 
given Initial data surface provided that certain conditions are met. Pre- 
vious investigators had used this approach in conjunction with an implicit, 
iterative finite-difference scheme to compute the supersonic flow over 
circular cones at angle of attack. In the present study, a new implicit, 
noniterative algorithm was developed which provides better computational 
efficiency than previous algorithms and is not restricted to conical 
shapes. 

The two approaches were used to calculate the supersonic flow over 
a cone at high angle of attack and a delta wing with a sharp subsonic 
leading edge. Both of the approaches accurately predicted the viscous 
and inviscld features of the two flowfields. In particular, the location 
of the main vortex above the delta wing was correctly predicted. 

During th-s study, an impoitant fact concerning the mathematical 
behavior ol the PNS equations wa3 discovered. It was found that the PNS 
equations remain hyperbolic-parabolic only if a fraction (u>) of the 
streamwlse pressure gradient is retained in the subsonic portions of the 
flowfield such as boundary layers. This fraction is given by 



1 + (y-l)M^ 


where is the local streamwlse Mach number. Using this knowledge, a 
technique was developed to prevent departure solutions which are attribut- 
able to the presence of the streamwise pressure gradient term in the 
streamwlse momentum equation. 


7 
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8. Development of a generalized PNS code 

The original parabolized Navicr-Stokes (PNS) computer code, described 

24 

in the previous section, was modified to permit the calculation of super- 
sonic viscous flows around arbitrary body shapes at high angles of attack. 
The previous code advanced the solution using axis-normal surfaces which 
allowed for most body geometries. However, for blunted bodies with large 
surface slope, the axial component of velocity in the inviscid part of 
the snock layer may be subsonic which prevented the previous code from 
computing this case. Therefore, it was necessary to generalize the code 
to march with solution surfaces that are more nearly normal to the local 
body surface. This was done in a perfectly general manner with a non- 
orthogonal three-dimensional coordinate frame. 

The generalized PNS code was used to compute the laminar flow over 
a 70° sweep slab delta wing having a spherical nose and a cylindrical 
leading edge. The body shape and flow conditions were chosen to match 
the experiments of Bertram and Everhart. Calculations were performed at 
angles of attack up to 41.5° and Mach numbers of 6.8 and 9.6. The com- 
puted shod, shapes, surface pressures, and heat transfer coefficients are 
in excellent agreement with the experiment. The results of these compu- 
tations are reported in Refs. 24-26. 

The complex turbulent flow in a three-dimensional corner was also 
computed using the generalized PNS code. Two simple eddy viscosity models 
were used for the turbulence modeling. The computed results were compared 
with the experiment of Peake and the complete Navier-Stokes computations 
of Horstman and Huns. Detailed comparisons of surface pressure, skin 




i / 
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friction, surface flow angle, velocity profiles and flow cn^le profiles 
show reasonable agreement. The results of these computations arc pre- 
sented in Refs. 27-28. 

9. Computation of the Space Shuttle Orblter flowfield 

The calculation of the complete inviscid-viscou3 flow around the 
Space Shuttle Orbitcr has been the aim of recent research efforts. While 
inviscld flow solutions and other solution techniques that couple the 
inviscid solution with the matching boundary-layer analysis have resulted 
in limited solutions for the windward side of Shuttle-like bodies, they 
have failed in general to predict the correct solution on the lee side 
and near the wing root region. In the present study, the laminar flow 
around the Space Shuttle Orbiter forebody was computed using the general- 
ized PNS code discussed previously. This generalized code is Ideally 
suited for the complex flow that develops around the leading edge and in 
the vicinity of the wing-body juncture of the Orbiter. The code will 
accept the initial solution on any starting data surface and it will 
allow the marching to proceed from one solution surface to another with- 
out any restrictions on the orientation of these solution surfaces. The 
only condition that must be satisfied is that the velocity component 
normal to the solution surface must be supersonic outside the boundary 
layer. 

The initial solution for the nose part of the Orbiter geometry was 
obtained with a 3-D time-dependent Navier-Stokes solver. It was necessary 
to employ a wind axis oriented coordinate system to obtain the initial 

29 

solution with the time-dependent code. A grid point clustering scheme 
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was employed to accurately describe the body shape by clustering points J 

i 

at the wing tip and at the wing-body juncture. The computed heat transfer i 

coefficients, pressure coefficients, and shock shapes compared favorably 
with the available experimental data for 0° and 30° angle of attack. The 

i 

results of these computations were reported in Refs. 26, 30, and 31. 

Work was initiated to modify the PNS code to permit real gas (equi- 
librium air) calculations. This modification was complicated by the 
presence of the Jacobians which must be evaluated for the finite-difference 
scheme. The real gas PNS code is currently being used to compute a Space 

Shuttle Orblter flowfield. The results of this computation will be coo- 

32 

pared with STS-3 flight data in a forthcoming paper. 

CONCLUDING REMARKS 

The research performed during this grant is described by the 32 
documents listed in the next section. This documentation consists of: 

8 Journal articles 

11 Technical papers presented and preprinted 
5 NASA Contractor Reports/Technical Memorandums 
4 Ph.D. dissertations 
4 M.S. theses 

In addition, the present grant supported the following students as 
Research Assistants at Iowa State University: 

. i 

1. Robert A. Mohling 

I 

2. Terry L. Holst : 

3. Paul H. Mugge 

% 

i 

\ 

I 

* 

* 


y 


* 


u 


J. Richard Eislcr 


5. Yvon C. Vigncron 

6. J. -Antoine Dcsideri 

7. Ethiraj Venkatapathy 

8. Iho K. Tree 

9. Dincsh K. Prabhu 

10. David C. Stookeaberry 

The late John V. Rakich of NASA Ames Research Center served as Technical 
Monitor for this grant until his death on April 9, 1983. The research in 
this grant was greatly enhanced by his inspiration and guidance. During 
the remainder of the grant period, Joe C. Marvin of NASA Ames Research 
Center served as the Technical Monitor. 
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